Neutron scattering in molecular liquids: 
Influence of orientational degrees of freedom and the prepeak in a fragile glass former 
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The intermediate scattering function S n (q,t) for neutron-scattering is expanded with respect to 
a complete set of correlation functions which describe the dynamical correlations in a molecular 
liquid. For the static ns-structure factor S n (q) of a system of diatomic molecules the results of 
the expansion are compared with the exact results from a MD-simulation and it is shown that the 
Sears-expansion, which is commonly used to interpret such data, fails in the supercooled regime. 
The representation for S n (?) is used to draw conclusions about the g-dependence and especially the 
origin of the prepeak. 

PACS numbers: 61.12.-q, 61.20.-p, 61.25.Em 



Neutron scattering is one of the most important tools 
to determine the structure and dynamics of condensed 
matter. For instance, it has been applied widely in the 
study of supercooled liquids and the glass-transition, (see 
e.g. ref. (lj and for more recent work [||^]), where it has 
been of great value to test theories such as the mode- 
coupling-theory (see e.g. HHI). One of the major ad- 
vantages of neutron scattering is that the neutron as an 
electrically neutral probe is not influenced by the electron 
cloud of the target atoms but interacts only with the nu- 
cleus. Consequently neutron scattering can be directly 
interpreted in terms of the atomic motion. However, if 
one wants to examine a molecular system not only with 
respect to the motion of the constituent atoms but in 
terms of the molecular units some care has to be taken. 
Since the superposition of scattering from atomic sites is 
in general not equivalent to scattering from the molecular 
center of mass, neutron scattering from molecular liquids 
is sensitive to orientational as well as translational cor- 
relations in the system. Attempts to account for this 
fact in the analysis of ns-data include the expansion into 
"partial waves" proposed by Sears ||. Originally con- 
ceived for liquids with negligible correlations between ro- 
tational and translational motion it was shown recently 
to describe well the incoherent scattering in supercooled 
water Jt[. The expansion proposed in the present Let- 
ter is closely related to that of Sears but does not make 
any assumptions on the strength of the coupling between 
rotation and translation. As a consequence it is able to 
describe well also the coherent scattering in the super- 
cooled regime. In addition it offers interesting new as- 
pects for the interpretation of the q-dependence of the 
static ns-structurefactor S n (q). 

A feature of special interest in supercooled liquids and 
glasses is the appearance of a prepeak in the static ns- 
structure factor at a g-value that corresponds to dis- 
tances larger than the average nearest-neighbour dis- 
tance. As a sign of intermediate range order prepeaks 
have been studied in a variety of systems (see e.g. || fuj). 
Mostly they have been attributed to the network struc- 
ture of strong liquids but as our present analysis of a 



system of rigid diatomic molecules shows prepeaks can 
appear even in the most simple molecular systems. 

We consider a set of N rigid molecules of identical ge- 
ometry each consisting of s atoms. The assumption of 
rigidity is one that is commonly used as well in the the- 
oretical analysis as in computer simulations of the liquid 
state. It is justified on the ground that at the tempera- 
ture of the liquid only the lowest vibrational states are 
populated. The starting point of our analysis is the fol- 
lowing site-site representation of the intermediate scat- 
tering function (cf. ||) 
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where a J h ( a i nc ) is the coherent (incoherent) scatter- 
ing length and Rj U (t) the position at time t of atom 
v in molecule j. The brackets (..) denote the canoni- 
cal average over initial conditions. Unlike the geometry 
which is identical for all molecules the scattering lengths 
are allowed to differ from molecule to molecule. This 
assumption is quite realistic since the chemical struc- 
ture is independent of the isotopic composition whereas 
the scattering lengths are. The translational and orien- 
tational motion are separated by introducing center of 
mass Rj(f) and relative coordinates Vj V (t), i.e. Rjv(i) = 
Rj(£) + *ju(t)- Analogously to Sears || we use the 
Rayleigh-expansion of the plane waves 



e iq-r,v(t) = ^ [4tt(2Z + 1)]* *' x 

In 



(2) 



to express the relative coordinates Yj V {t) by the polar co- 
ordinates (r u ,6 v ,4> u ) of atom v with respect to a body 
fixed frame with origin at the center of mass of molecule 
j and by the Euler-angles flj{t) denoting the orienta- 
tion of the body-fixed system with respect to the labo- 
ratory frame. 
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Without loss of generality the z-axis of the laboratory 
system has been chosen to point in direction of q. The 
special functions appearing in eq.(^) are the spherical 
Bessel functions ji, the spherical harmonics Y" 1 and the 
Wigner functions T> l mn . Conventions are chosen accord- 
ing to the textbook of Gray and Gubbins Oj. The 
range of indices is 1,1' = 0,1,2,...; — I < m,n < I; 
—I' < m, n' < I'. Thus we establish a connection between 
the intermediate scattering function S n {q,t) of neutron 
scattering and the set of functions 



S^Aq,t) = i 1 '- 1 W + l)(2l> + 
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which can be split up into self (j 
(j ¥= f) giving 

->m(s) 



j') and distinct part 
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This set of correlation functions is the generalization to 
arbitrary molecules of the correlators recently introduced 
for a single linear molecule in an isotropic liquid |l5| ] and 
for molecular liquids of linear molecules |l4|]. Similar 
correlation functions have been used previously in the 
study of molecular liquids |l^,[l5| . The tensors @ have 
a number of symmetry properties which are discussed in 
refs. [ |l6|]l4f| . We just want to point out that S® 00 (q, t) is 
the usual intermediate scattering function for the center 
of mass whereas 5 Z ™ i/ n ,(q, t) for I or I 1 different from one 
reflect correlations between higher "multipoles" of the 
microscopic density. 

Inserting eq.(||) into ([[]) one obtains with (Q) the final 
result: 
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with the coefficients 



(6) 



(7) 
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Here x = l/N^j x 3 denotes the average over molecules. 
A similar result for diluted linear molecules in an 
isotropic liquid is given in ref. Ej ]. In contrast to the 
original approach by Sears we will not invoke factor- 
ization of correlations between (i) all translational and 
rotational motion and (ii) the orientational degrees of 



freedom of different molecules. This asumption which 
is quite common in the analysis of neutron as well as 
light-scattering turns out to be not quite satisfactory in 
the supercooled regime as we will demonstrate in the 
following. The expansion (Q) has a number of interest- 
ing properties. Immediately obvious is that eq.(^) can 
be Fourier-transformed to give a corresponding relation 
between the spectra or susceptibilities Further, we no- 
tice that only correlators with m = enter into S n (q 7 t) 
(which is a consequence of the isotropy of the fluid) and 
that the coefficients bi n ^ n i (q) are completely determined 
by the geometry of the molecule via (r u ,9 u , (f) u ) and av- 
erages of the scattering lenght of the constituent atoms. 
The dependence on the wavevector q enters through the 
spherical Bessel functions ji(qr v ) and is thus connected 
to the radii r u and the value of I. This connection be- 
tween q and I will be very useful to justify the cutoff 
of the summations over I in eq. (||) and to analyse the 
ns-structure factor. 
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FIG. 1. Static molecular correlation functions up to 
1,1' < 2 for the lowest temperature T = 0.48 (in units of 
the Lennard-Jones energy of the A-atoms) of the simulation 

To demonstrate the quality of expansion (|^) we con- 
sider the special case of diatomic rigid molecules. From 
a recent MD-simulation [0 we have taken the static 
molecular correlation functions according to eq.(|§|) and 
we have evaluated the exact result for the static ns- 
structurefactor using eq. (0) . In that simulation a system 
of 500 molecules was considered each consisting of two 
atoms labeled A and B. The molecular interaction was 
represented as the superposition of Lennard-Jones po- 
tentials between the atoms. The molecular bond length 
was fixed at d = 0.5 in units of the Lennard-Jones radius 
of the A-atoms which we will use throughout this Let- 
ter. Further details about the simulation can be found in 
references JTt| . Since for linear molecules the third Euler- 
angle \ doesn't play a role the distinct parts of the molec- 
ular correlation functions are nonzero for n = n' = only 
and the summations over n, n' can be carried out for the 
self part. Therefore we will skip the indices n, n' as well 
as the superscript m — in the following. The static cor- 
relation functions Su>(q) evaluated up to 1,1' < l co = 2 
are shown in Figure |l| for the lowest temperature of the 
simulation. For > l co only the self part sff (q) = 5u> 
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is considered in the following calculations. Specializing 
to linear diatomic molecules and the static case eqn.(|[|^) 
simplify to: 

Sn(q)= £ b% h (q)S$(q)+b inc ( q ) (8) 

l,l'<l c o 

6^( 9 ) = [(2/ + l)(2/' + l)] i i t (y)jV(y) x (9) 

X \ (<~H+ (-1)^) (<H + (-1)^) 

& mc w =^(^+^+^ii+<:)+ (io) 

a coh a cohio{qd) 

The reader should notice that l co = just yields the 
usual Sears-expansion with the molecular structure fac- 
tor b mc (q) and the intcrmolccular form factor b'^ l (q) for 
uncorrelated rotation. In Figure || we show the compar- 
ison between the exact S n (q) and the expansion (||) for 
lco = (Sears), l co = 1 and 2 at the lowest tempera- 
ture of the simulation. The expansion for l co = 2 is in 
perfect agreement with the exact S n (q) up to q « 9 and 
in qualitative agreement up to almost q w 18 whereas 
the Sears-expansion shows hardly more than qualitative 
agreement up to q « 8. If the scattering length of A- and 
B-atoms are equal it follows from (^|) that only even I and 
V contribute to S n (q). Hence S n (q) from (|8|) with l co = 1 
is identical to the Sears-result. If neutron-scattering can 
distinguish between A- and B-atoms the l co — 1-result 
is still not much better than the Sears-expansion (cf. 
Fig. (Q)) except for the prepeak at q w 3 which does not 
show up in the latter (see inset of Fig. (|2|)). This failure 
of the Sears-expansion clearly shows that in the super- 
cooled regime the assumption of uncorrelated rotation of 
different molecules does not hold. Upon increasing the 
temperature these correlations become weaker and we 
find that the contributions coming from Sw(q) with 
different from get less important. Whereas usual ways 
to account for orientational correlations are restricted to 
assuming specific orientational models like in the analysis 
of Dore et al. (l8) or to fitting the data under the assump- 
tion of a preferred orientation |]l9f our approach is an ex- 
act expansion (provided l co — oo) that makes no ad hoc 
assumptions about orientational correlations. Of course 
the calculation of Vn , (q) can not be done exactly, but 
requires either approximation schemes [|l2| , p^[ or simula- 
tional methods. In the case of the time-dependent corre- 
lators 5*;™ ;/„/(<?, t) one may use the results following from 
the molecular mode coupling theory |l^] . 

Comparing the q-dependence of the input-data in Fig- 
ure [l] with that of S n (q) reveals the origin of the different 
peaks in the ns-data. Whereas the position of the pre- 
peak coincides with the first peak in Sn(q) the main peak 
has its origin in the center of mass correlations Soo(q) 
though it is slightly shifted to higher q due to the con- 
tributions from So2(q) and 522(g)- The peaks at q w 13 



and q ~ 18 are easily recognized as stemming mainly 
from the corresponding ones in the correlator S , 22(<z). 



S n (q) is 
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FIG. 2. Comparison of the different approximations to 
S n (q) at T = 0.48 for at h = 1.4, a B coh = 0.25, af nc = af nc = 0. 
The inset shows an enlargment of the prepeak region. 

These obvious identifications are supported even more 
by a closer inspection of the coefficients 6^°' l (g) which 
will also explain why the structure in S n (q) at q > 18 is 
missed by the cut-off we have chosen. As stated in the 
discussion of the general result (||) the g-dependence of 
bw(q) is closely related to the values of I and /' through 
the Bessel functions These are oscillating func- 

tions that decay to zero for increasing argument pOj . 
Whereas jo(x) has a maximum at x = the functions 
jl (x) with / > 1 start at zero and the position of the first 
maximum shifts to higher x with increasing 1. This prop- 
erty is inherited by the coefficients bf°, h (q) and leads to 
the consequence that correlators with small I give a con- 
tribution at low q whereas big I determine the behaviour 
at high q (apart from the incoherent contribution from 
b mc (q)). This is also the reason why the expansion with 
l co = 2 works less good for q > 9 and particularly for 
q > 18 since higher values of I and I' would be needed. 
A quantitative analysis of the coefficients at the position 
q m 3 of the prepeak shows that about 80% of the prepeak 
amplitude originates from Sn(q). Therefore the appear- 
ance of the prepeak will depend sensitively on the ratio 
of the scattering lengths, as demonstrated in Figure |[ 

In the final part of this Letter we want to analyse fur- 
ther the appearance of a prepeak in S n (q) for diatomic 
LJ-molecules. As we have already pointed out in the dis- 
cussion of the g-dependence of S n (q) the prepeak is con- 
nected with the main maximum of the correlator Sn(q). 
The center of mass correlations show no structure at all 
at q ~ 3 which is a remarkable contrast to the situation 
encountered in supercooled water where the intermediate 
range order shows up as a prepeak in Soo(q) pi] - This 
indicates that different mechanisms are responsible for 
its existence in the present fragile glassformer and the 
network-former. Taking into account the shape of the 
spherical harmonic function Y® oc Dg one can conclude 
that Sii(g) can be viewed as the correlations between 
molecules with parallel orientation (of next nearest neigh- 
bours). Thus the connection of intermediate range order 
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with Sxx(q) is in good agreement with the findings of Mi- 
sawa JjlJ who proposed a staggered parallel orientation 
as the preferred orientation of nearest neighbours in liq- 
uid halogenes. In another paper Misawa j|] showed that 
the existence of a preferred orientation can give rise to 
a prepeak depending on the ratio of interatomic spacing 
and bond length. 



S„(q) 0.2 



FIG. 3. Dependence of the prepeak amplitude on the ratio 
of coherent scattering lengths of A- and B-atoms. 

In conclusion we can say that the proposed expan- 
sion for the intermediate scattering function of neutron- 
scattering offers a general way of describing orientational 
correlation effects. It indicates interesting connections 
between the g-dependence of S n (q), the orientational cor- 
relations for different I and the molecular geometry which 
are useful in the interpretation of neutron-scattering re- 
sults for molecular liquids. In particular, we have shown 
that orientational correlations play a crucial role in the 
supercooled regime. This is consistent with the recent 
results by Bermejo et al. ]22| . These correlations are 
also responsible for the existence of a prepeak, whose in- 
tensity depends both on the scattering lengths and the 
temperature. 
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